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Abstract 

Short distance static quark-antiquark interaction is studied sys- 



o 
m 

co 

>. 

^-C) , tematically using the background perturbation theory with nonper- 

turbative background described by field correlators. A universal linear 
^vi ■ term c ^ ar is observed at small distance r due to the interference be- 

C tween perturbative and nonperturbative contributions. Possible mod- 

O^ ' ifications of this term due to additional subleading terms are discussed 

On ' 

--^ ■ and implications for systematic corrections to OPE are formulated. 

ff- 1 Introduction 



It is more than 20 years ago that the power correction has been computed in 
OPE [1] laying ground for numerous later applications in QCD. Since then 
OPE is the basic formalism for study of short-distance phenomena, such as 
DIS, e + e~ annihilation and, with some modifications, heavy quark systems. 

Interaction of static charges at small distances has drawn a lot of atten- 
tion recently [2-4]. The theoretical reason is that the appearance of linear 
terms in the static potential V(r) = const r, where r is the distance between 
charges, implies violation of OPE, since const ~ (mass) 2 and this dimension 
is not available in terms of field operators. There are however some analytic 
[5,6] and numerical arguments [7,8] for the possible existence of such terms 
0{m 2 /Q 2 ) in asymptotic expansion at large Q. 

Of special importance is the sign of the mass squared term. It was argued 
recently in [3] that the small distance region may produce tachyonic mass 
correction and this correction was studied selfconsistently in different QCD 



processes. In particular the correct (positive) sign of linear potential at small 
distances comes from tachyonic gluon mass, while positive mass squared term 
produces negative slope of linear potential. 

On a more phenomenological side the presence of linear term at small 
distances, r < T g , where T g is the gluonic correlation length [9,10], is required 
by at least two sets of data. 

First, the detailed lattice data [11] do not support much weaker quadratic 
behaviour of V(r) ~ const r 2 , following from OPE and field correlator method 
[9,10], and instead prefer the same linear form V(r) = or at all distances (in 
addition to perturbative — ^°s. term). Second, the small-distance linear 
term is necessary for the description of the fine structure splittings in heavy 
quarkonia, since the spin-orbit Thomas term V t = -^^ is sensitive to 
the small r region and additional linear contribution at r < T g is needed 
to fit the experimental splittings [12]. Moreover the lattice calculations [13] 
display 1/r behaviour for the spin-orbit potential V[ in all measured region 
up to r = 0.1/m. 

In what follows we display the basic dynamics which produces tachyonic 
gluon mass and estimate its magnitude. 



2 Background perturbative theory 

In this letter we report the first application of the systematic background 
perturbation theory [14] to the problem in question. One starts with the 
decomposition of the full gluon vector potential A^ into nonperturbative 
(NP) background B^ and perturbative field a M , 

A^ = B^ + a M , (1) 

and the 'tHooft identity for the partition function 

Z = J DA,e' s ^ = 1 J DBXB) J Da,e~ s ^ (2) 

where rj(B) is the weight for nonperturbative fields, defining the vacuum 
averages, e.g. 

< F*„(x)® B (x, y)F*,(y) > B = ^{5, x 5 va - 5,J uX )D(x - y) + A, (3) 



where F^ u , $ B are field strength and parallel transporter made of B^ only; 
Ax is the full derivative term [9] not contributing to string tension cr, which 
is 

a = J- f d 2 xD(x) + 0(< FFFF >) (4) 

The background perturbation theory is an expansion of the last integral in 
(2) in powers of ga^ and averaging over B^ with the weight rj(Bfj), as shown 
in (3). Referring the reader to [14] for explicit formalism and renormalization, 
we concentrate below on the static interquark interaction at small r. To this 
end we consider the Wilson loop of size r xT, where T is large, T — ► oo, and 
define 

< W > B ,a=< Pexpig / (B^ + a^dz^ > B ,a= exp{-V(r)T} (5) 

Expanding (5) in powers of ga^, one obtains 

<W >=W + W 2 + ...- V = V (r) + V 2 {r)+V 4 {r)+, (6) 

where V n (r) corresponds to {ga^) n and can be expressed through D,Ai and 
higher correlators [10,14] and its behaviour at small r is [9] 

V (r)=C r 2 + Cy + ..., r<T g ,V (r)=ar, r > T g (7) 

where coefficients are integrals of field correlators over Euclidean time, 

C = 0(< FF>)= / D(v)dv + O(Ai), C' = 0(< F 4 >). 
Jo 

It is this small r behaviour (7) which causes phenomeno logical problems 
mentioned above. 

Coming now to V^fV), describing one exchange of perturbative gluon in 
the background, one finds from the quadratic in a^ term in S(B + a) in the 
background Feynman gauge the gluon Green's function 

G lu , = -(I%8 lu , + 2igF*)- 1 , DT = d x 5 ca + gf cba B b x (8) 



The term W 2 can be written through G^ u as 

W 2 = g 2 dx 4 / dy 4 < Pexpiig I Badzu)Gu{x,y) > B (9) 

Jo Jo Jc 



where for G44 one can use the Feynman-Schwinger representation (FSR) [10]. 
The simplest form of FSR obtains when one can neglect or expand in powers 
of gluon spin interaction (paramagnetic term 2igF B in (8)). 
Doing this expansion, G^ can be written as 

G = -D~ 2 + D- 2 2igF B D- 2 - D- 2 2igF B D~ 2 2igF B ' D~ 2 (10) 

the first term on the r.h.s. of (10) corresponds to the spinless gluon exchange 
in the background B^, which can be written using FSR [10] as 

G u (x : y)= ds(Dz) xy e~ K Pexpig B^dz^, K = -\ z 2 dr (11) 

jo Jc(z) 4 Jo 

Here B^ is in the adjoint representation of SU(iV c ). At large N c one can 
use the 'tHooft rule to replace the gluon adjoint trajectory C(z) by the double 
fundamental trajectory which forms together with the original rectangular 
contour C in (9) two closed Wilson loops (see [14] for details and discussion). 

The average over B^ in (9) then reduces (again in the large N c limit) to 
the product of two averaged Wilson loops, namely 

W 2 = / ds(Dz) xv e- K dx A dyA < W(d) >< W(C 2 ) >, K=- z 2 dr, 
Jo 4 Jo 

(12) 

where C\ and C 2 are two contours obtained from the rectangular Wilson 

loop when two points on it, x and y, are connected by a double line of 

gluon trajectory. It is convenient to choose the surfaces Si inside C\ and C 2 

as consisting of two adjacent pieces Si{Ci) = Si{plane) + S(A), one lying 

on the plane of original Wilson loop and another piece A, perpendicular 

to the plane and bounded by the trajectory. Using now the nonabelian 

Stokes theorem [9] and cluster expansion for the average < W(Cj) >, one 

obtains that bilocal correlator of fields with points on the two pieces vanishes 

since < EiEk >~ Sik, and < EiBk >= 0(Ai) and vanishes by symmetry 

arguments. Trilocal and higher correlators for dimensional resons bring with 

them higher powers of distance r and can be neglected. 

Hence one has 

< W(d) >< W(C 2 ) >=< W A > 2 W (13) 

where 

W A = ex P(-7jjj- J D ( x ~ y)d^nv{x)da^ v {y)) (14) 



Two different regimes are possible for (14). In the small distance region, 
r < T g , the sizes of the surface A are of the order of r and one can replace 
D(x — y) — > -D(O) in (14), and for dimensional reasons the only possible 
contribution is 

W A = l + O(D(0)r 4 ) (15) 

Thus one obtains a correction 0(r 4 ) to the perturbative potential 1/r, and 
hence no linear term. 

In the large distance region, r ^> T g , one obtains the area law for Wa, 

W A =* exp(-aS A ) (16) 

Insertion of (16) in (12) yields a massive propagator of a spinless hybrid with 
mass m = m A at large r, which corresponds to the first excitation of the 
open string with fixed ends. 

Summaring one can rewrite W 2 as 

W 2 = W J dx 4 J dy 4 < G(x, y) > (17) 

where G(x,y) > is the Green's function of the spinless hybrid. One can 
satisfy the properties (15), (16) representing < G > as the propagator of a 
particle with variable mass mo(p), 

< G (P) >= 2 j 2? N ( 18 ) 

where ml(p) < 0(\),p — >■ 00, and ml(p — >• 0) = m A . 

Thus at small distances (large p) < G(p) > describes the usual massless 
gluon exchange, whereas at large distances it describes the propagation of 
the spinless hybrid. 

Hence at small r < T g the background field B^ in D~ 2 is not operative 
and one can replace D~ 2 by the free gluon propagator d~ 2 . The contribution 
of the second term to W2, dx tl D~' 2 F^ u D~ 2 dy u , obtains when x and y are on 
adjacent sides of the Wilson loop and therefore does not affect V(r). In what 
follows we concentrate on the third term in (10), W 2 , 

W 2 (3) = 4g 2 ( dx 4 f dy 4 < G(x,u) > d\ < g 2 F£(u)F%(v) > 
Jo Jo 

x < G(u,v) > d 4 v < G(v,y) > (19) 



One can rewrite (19) as 

WP = 2g 2 / dx, / dy, / «J?L———\*(k\ (20) 



Jo 4 7o y *J (P + mg(fc)) 2 (27r) 4 
where we have defined, having in mind (4) 

2/IX r D(z)e~ ikz d 4 z ,, , 3aN c . , 

"< t) = 3 / w ^ (0 > = — ■ (21 » 

Doing integrals over dx^dy^, one gets 

r(3) _ T f d 3 ke l ^a s {k)X 2 {k) 



W * = ^ / 7t2 I 2/, ^2 = " A ^ r ) T (22) 

?r J (k 2 + m 1 (k)) 2 

To estimate the integral (22) one can take a s (k)X 2 (k) out of integral at some 
effective point k and calculate the rest in a simple way, assuming m to be 
constant, and expanding result at small r. 
In this way one obtains 

AV 2 (r) = ^M^MJ_ e -m r = as{ k )^k )(-±- + r + 0(r 2 )) (23) 
rrior orriQ m® 

Analysis of (21) tells that A 2 (/c) is a rather weak function of the argument, 
and to get an idea of the magnitu 
A 2 (0) = ^ yielding for AV 2 (r) 



and to get an idea of the magnitude of AV 2 (r), one can approximate A 2 (A;o) — 

3crN c 



AV 2 (r) , r < T g . (24) 

Note, that had we renormalized a s in (23), (24) in the standard way we 
would meet the IR divergence of the running a s (k), since the corresponding 
momentum k$ is in the IR regime for the constant term ( — — ) in (23). 
However, the NP background formalism predicts IR modification of a s (see 
[14,15] for details and discussion), the so-called freezing a s behaviour, which 
from heavy quarkonia fitting was found in [16] to yield maximal a s at small 

k, 

a s (max) = 0.5 (25) 

For the linear term in (23) the situation is different^ and the effective value 
of k is of the order of 1/r. Hence it is more appropriate to present (23) in 
the form 

AF 2 (r) = a s (l/r)/i 2 (l/r)r (26) 

1 The author is grateful to V.I.Zakharov for the discussion of this point 



One may wonder whether behaviour (24) holds also at large r, thus in- 
creasing linear potential Vo(r). However the form (24) was obtained at small 
r < T g , while at large r next terms of expansion in (10) are contributing and 
the whole series (10) should be summed up explicitly. 

One can perform the summation, replacing first for simplicity D — »■ d 
in (8), (small r or relatively large q are considered) and assuming that only 
bilocal correlators of F^ u are nonzero. One obtains 

G = -(d 2 + 2igF)- 1 ^G(k) 



k 2 -^{kY 

where ^ 2 {k) > 0. Thus G(k) acquires a pole at real value fi 2 (k) in Euclidean 
space-time, signalling appearance of a tachyon. From physical point of view 
this result is a consequence of paramagnetic attractive interaction of gluon 
spin with NP background, yielding negative correction to the gluon selfen- 
ergy. A similar term occurs for a quark due to its spin interaction with 
background [17]. 

Recently a negative (tachyonic) mass shift was observed due to the gluon 
interaction with the stochastic background in [18]. 

It is meaningful, that the same paramagnetic term F^ v in G^ u yields 
negative contribution to the charge renormalization (asymptotic freedom) 
[19]. In fact negative paramagnetic effective action S^J [19] gives rise to 
the negative (tachyonic) mass since both are connected, — n 2 5 i ' " 



Therefore one may expect that the phenomenon of tachyonic gluon mass is 
pertinent to nonabelian theories. 

The existence of the tachyon reveals paramagnetic instability of the ob- 
ject (gluon or quark), if stabilizing mass is not created by some additional 
mechanism. In our case this mechanism is the creation of hybrid mass due 
to the same confining correlator D(x) when D 2 is used and not d 2 in (8). As 
was explained above, the hybrid mass is created at larger distances, r ^> T g , 
so that for illustrative purposes (referring the hybrid mass to the gluon in 
question), one may write the total gluon propagator G(k) as 

G(k) 



k 2 + m 2 {k) 



where m 2 (k) = m^fc) — n 2 (k), and rn^{k) is dominating at small k (large 
distances) while n 2 (k) dominates at large k (small distances). A simple 
example is provided by m 2 {k) = /x 2 ^ 2 ~ 2 , in which case one can calculate 



/* 



gluon exchange potential V(r) = — J -4^e lkr G(k) explicitly to yield at small 

r: 

1 u 2 u 2 

V(r)~-- + ^-r, AV-^-r 

in agreement with the form (24). At large distances one should calculate the 
exact gluon Green's function (8) in the Wilson loop W^ 2 \ which describes to 
the propagation of the hybrid state with two static quarks QQ at the ends 
of the string. 

Such a state was considered both analytically [20] and on the lattice [21] 
yielding excitation energy (which corresponds to the mass m(0)) around 1 
GeV. 

Hence (24) is only a small distance approximation of the hybrid exchange 
potential, where the dominant paramagnetic contribution is kept in the ef- 
fective gluon mass. 



3 Other possible corrections to V(r 

In doing perturbative expansion in (2) one encounters other terms in S(B+a) 
which potentially yield interference contributions of perturbative a p and NP 
B^ fields. Of special importance is the term L 1; 

L x = j a v D,{B)F B u d A x (27) 

Correction due to (27) in the gluon propagator was studied in Appendix 1 of 
[14] and can be written in the form (20) where fi 2 {k 2 ) — »■ J\ is now expressed 
as 

J x (k 2 ) = fd 4 ze~ zkz < DF B (z)DF B (0) > (28) 

Now the integrand in (28) can be written as a sum of two terms [22], 

< D p F*(z)D x F^(u) >= d p d x < F B (z)F B (u) > +0(< FFF >) (29) 

The first term on the r.h.s. of (29) yields 5J\ ~ k 2 and hence only a NP 
correction to the Coulomb term, while the second term, < FFF >, would 
give, using dimensional arguments, correction AV^ ~< FFF > r 5 , negligible 
at small r. 



Another type of correction occurs from the interplay of multiple color 
Coulomb exchanges and one NP contribution, considered in [23]. This effect 
can be accounted for by the replacement 

Not 
D(x, t) -► D(x, t)exp{ pi) (30) 

Here t is the Euclidean time and the exponent in (30) accounts for the 
difference of potential in singlet and octet channels. Insertion of (30) in (7) 
yields an additional suppression of the r 2 dependence. This result coincides 
with the correction obtained in [24] in a different way. 

Finally we consider in this section the correction due to the freezing be- 
haviour of the coupling a s [14,15], 

a f(r) = *L - a (0)( r ) _ ai " )m B r2 m 

As was noticed by F.J.Yndurain [2], expansion of the freezing Coulomb 
potential yields a linear term 



C 2 a(°) 



mlr 



AV c (r) = + ~*,'jy +■■■ (32) 

Here tub is the double hybrid mass, (tub = l.lGeV from the fits to experi- 
ment [16]), and the term (32) is always much smaller than the Coulomb term 
because of condition m B r 2 <ti 1. One may consider AV C in (32) as coming 
from the additional 1/p 2 in a s (p) , as was suggested in [8], but here there is 
the log term in the denominator, reminding that the pole is coming from the 
expansion of the freezing a s (p). 

One should note that there is no double counting in adding (24) and (32), 
since (24) is obtained from the one-gluon exchange (OGE) process, while (32) 
is due to the one-loop corrections to OGE. However the region of validity of 
(32) is always smaller than (24). 

4 Discussion and conclusion 

The analysis done heretofore concerns static interquark potential and re- 
veals that even at small distances NP background ensures some contribu- 
tions which come from relatively small intermediate distances, I < T g , and 

9 



encoded in the mass squared term 

H 2 = JL-l, n » 0.5GeV (33) 

which generates the potential (24). 

Applying the same NP background formalism to other processes of inter- 

2 

est, one would get similar corrections of the order of ^ , as e.g. in OPE for 
e + e~ annihilation. 

One example of this kind is the calculation of the field correlator 

< F liV (x)^(x, y)F\ a (y) >. The leading contribution can be written using the 
gluon propagator G (8), 

< F$F >= d^d x G ua + perm + G^ x G ua + perm (34) 

where we have suppressed $ and perm, denotes terms obtained by the permu- 
tation of indices with the proper change of sign. Insertion of expansion (10) 
into (34) yields in addition to the standard perturbative term 0((x — y)~ 4 ) 
a contribution proportional to (33), 

< F(x)$F(0) >~ ^ + ^ + ... (35) 

where both C\ and C2 are positive computable numbers. 

A recent lattice study of a similar quantity [7] as a function of UV cut-off 
A reveals the possible presence of the 0(A 2 ~ x~ 2 ) term. 

It is clear that appearance of /x 2 , which is an integral of nonlocal entity 

< F(x)F(0) > over a NP scale, x ~ T g , violates the original OPE of Wilson 
[25], proved in the pure perturbation theory, and the extended OPE of Shif- 
man, Vainshtein and Zakharov [26], where NP contributions enter as matrix 
elements of local operators. 

This extended form of OPE can be considered as a physically motivated 
assumption, and an explicit treatment done here within the background per- 
turbation theory (BPT) reveals that some extra terms should be added to 
OPE, the first of which, fi 2 /p 2 , was discussed in [3]. 

One might ask at this point, how rigorous and selfconsistent is BPT, 
with nonperturbative background given by correlators. One should stress 
here, that BPT is a consistent and systematic method, but not a complete 
one, since no recepee was suggested above for calculation of NP correlators, 
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and the NP configurations are introduced by hand just as it is done in QCD 
sum rules [26]. 

However recently the situation has changed. In [27] equations have been 
derived in the limit of large N c for vacuum correlators, from which correlators 
can be computed one by one explicitly. In this way the exponential form 
of the lowest correlator, < F(x)<&F(0) > was defined analytically [27] in 
agreement with lattice studies, and analytic connection between T g and o 
was found, yielding T g in a good agreement with lattice data [28] . 

The main conclusion from these studies is that NP configurations ap- 
pear as a selfconsistent solution of nonlinear equations which violates spon- 
taneously scale symmetry pertinent to these equations and their perturbative 
solutions. 

/,From this point of view the NP background exploited here can be iden- 
tified with scale violating NP solutions in [27], and the BPT method is made 
complete. 

On the phenomenological side the account of the correction (24) and (32) 
in the total potential 

V(r) = V (r) - ^M + AV 2 (r) 
r 

may provide for V(r) a simple "linear plus Coulomb" picture for all distances 
which is in better agreement both with experiment [12,17] and with lattice 
data [11, 13]. 

To conclude: in this study perturbative - nonperturbative interference 
was shown to provide additional OPE terms, absent in the usual local OPE 
form. In addition there are purely NP contributions [14] which are also 
outside of the standard lore, and will be discussed elsewhere. 

The author is grateful to V.A.Novikov and V.I.Shevchenko for fruitful 
discussions, and to V.I.Zakharov for discussions, correspondence and very 
useful remarks. 

The financial support of RFFI through the grants 97-02-16406 and 97- 
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